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We demonstrate that the Higgs mechanism in three-dimensional topological superconductors exhibits unique 
features with experimentally observable consequences. The Higgs model we discuss has two superconducting 
components and an axion-like magnetoelectric term with the phase difference of the superconducting order 
parameters playing the role of the axion field. Due to this additional term, quantum electromagnetic and phase 
fluctuations lead to a robust topologically non-trivial state that holds also in the presence of interactions. In 
this sense, we show that the renormalization flow of the topologically nontrivial phase cannot be continuously 
deformed into a topologically non-trivial one. One consequence of our analysis of quantum critical fluctuations, 
is the possibility of having a first-order phase transition in the bulk and a second-order phase transition on the 
surface. We also explore another consequence of the axionic Higgs electrodynamics, namely, the anomalous 
Hall effect. In the low frequency London regime an anomalous Hall effect is induced in the presence of an 
applied electric field parallel to the surface. This anomalous Hall current is induced by a Lorentz-like force 
arising from the axion term, and it involves the relative superfluid velocity of the superconducting components. 

The anomalous Hall current has a negative sign, a situation reminiscent of, but quite distinct in physical origin 
from the anomalous Hall effect observed in high-superconductors. In contrast to the latter, the anomalous 
Hall effect in topological superconductors is non-dissipative and occurs in the absence of vortices. 

PACS numbers: 74.20.-z,74.25.N,03.65.Vf,11.15.Yc 


I. INTRODUCTION 

Topological solid states of matter 1 2 have the property that 
their bulk states are gapped, while states at the boundaries are 
gapless and protected by some discrete quantum symmetry. 
The topological aspect emerges when considering the trans¬ 
port properties of the boundary states, where the transport 
current happens to also be a topological current. The most 
well established topological solid states of matter are topolog¬ 
ical insulators (TIs), which are gapped in the bulk and have 
helical (or chiral) gapless states at the boundaries which are 
protected by time-reversal symmetry. Helical here means that 
the electronic spin is locked to momentum due to strong spin- 
orbit coupling. Thus, the boundary states have an helicity 
determined by the eigenvalues of cr • k/k at each boundary. 
Topological insulators have been predicted to exist 3 and con¬ 
firmed experimentally in subsequent papers 4 . Although many 
of the materials investigated experimentally are not perfect in¬ 
sulators in the bulk, the observed boundary helical states are 
robust features of these materials. 

There exists another predicted topological solid state of 
matter, namely topological superconductors (TSCs) 1 2 , where 
the experimental situation is less clear. TSCs follow a 
symmetry classification scheme closely related to TIs, as 
far as Hamiltonians of the Bogoliubov-de Gennes type are 
concerned.^ Just like TIs, TSCs have gapped states in the bulk 
and symmetry-protected gapless states at the boundaries. Un¬ 
like TIs, in TSCs the 7/(1) symmetry is broken, either spon¬ 
taneously or by proximity effect. The gapless boundary states 
are Majorana fermions, which are fermionic particles that are 
their own anti-particles. In order to support such states at the 


boundaries, the topological superconductivity must feature a 
p-wave type of gap. Particle-hole is the underlying symmetry 
protecting the boundary Majorana states. In one dimension, 
a paradigmatic simple model for topological superconductiv¬ 
ity has been proposed by Kitaev 7 where the Majorana zero- 
energy states live at the ends of a quantum wire. An experi¬ 
mental way of realizing a superconducting state in a quantum 
wire is by proximity effect. In this case a semiconducting wire 
with strong spin-orbit coupling is deposited on the surface of 
an s-wave superconductor in the presence of an external per¬ 
pendicular magnetic field. Then by proximity effect p -wave 
like superconductivity is induced on the wire for a certain 
range of parameters P - 9 There are some experimental signa¬ 
tures of Majorana modes in Indium antimonide nano wires in 
contact with normal and superconducting electrodes.® More 
recently, strong support for Majorana boundary zero modes 
have been reported in an experiment with a ferromagnetic 
chain of Iron fabricated on the surface of lead. 11 

Three-dimensional (3D) TSCs have also been discussed 
theoretically, in particular focusing on vortex physics 12 and 
possible topological phase transitions! 12 A distinctive feature 
of both three-dimensional TIs and TSCs with respect to their 
non-topological counterparts is the topological magnetoelec¬ 
tric response induced by a mechanism similar to the so-called 
chiral anomaly. 13 When fermions in topological materials in¬ 
teract with the electromagnetic field, a Berry phase mixing 
electric and magnetic fields is induced.®® In TIs this occurs 
due to strong spin-orbit coupling that locks spin to momen¬ 
tum. The resulting Berry phases combine in the form of a 
so-called axion term, which is a magnetoelectric term ~ E • B 
with a periodic field, 6 , appearing as a coefficient.® This co- 


2 


efficient corresponds to a topological invariant implied by the 
Chern number of the band structure. If TR symmetry is pre¬ 
served, there are only two possible values for 6 in a 3D in¬ 
sulator, namely, 6 = 0 and n, the former value correspond¬ 
ing to a topologically trivial insulator.® In the case of TSCs a 
topological magnetoelectric contribution also arises, but now 
6 corresponds to the phase difference between order parame¬ 
ters of opposite chirality. 15 Recently, axionic superconductiv¬ 
ity has been also discussed in the context of doped narrow-gap 
semiconductors 

The 3D TSC constructed by Qi et al. 1 ^ should actually cor¬ 
respond to an interacting topological state of matter, going 
beyond the the symmetry classification of topological non¬ 
interacting theories for insulators and superconductors P Ac¬ 
cording to the standard symmetry classification and an argu¬ 
ment following from the gravitational anomaly,® 3D TSCs 
having TR invariance are classified by an integer topologi¬ 
cal invariant, belonging to the Dill class in the free fermion 
classification table. 6 Thus, unlike TR invariant TIs, the elec¬ 
tromagnetic axionic response of 3D TSCs would involve an 
axion field having values 6 = 0 or n, mod 2n. However, 
as pointed out in Ref. 15, the axion in TSCs is a dynami¬ 
cal phase variable associated with the superconducting order 
parameter, or Higgs field in the language of quantum field the¬ 
ory. Furthermore, the Meissner effect gives a mass to the pho¬ 
ton via the Higgs mechanism. This is an important difference 
with respect to the axion electrodynamics of TIs, where no 
U(l) symmetry is broken and the gauge field remains gap¬ 
less. Furthermore, nontrivial Chern numbers are associated 
with the different phases. For the simpler case featuring two 
Weyl fermions with opposite helicity, we have opposite Chern 
numbers, leading to a topological invariant given by the sign 
of the gap amplitude. 15 Thus, the arguments in Ref. 15 seems 
to point out to a Z 2 classification. However, since the rele¬ 
vant symmetry in the problem is U(l) x Z 2 , it has been shown 
recentl)®^that in topological superconductors with TR sym¬ 
metry the Z classification is reduced to Z\§. 

In this paper, we investigate the Higgs mechanism and 
anomalous Hall effect of three-dimensional TSCs within the 
model introduced recently in Ref. [15] In the simplest case, 
the model features two superconducting order parameters as¬ 
sociated with left and right fermion chiralities interacting with 
the electromagnetic field, and a topological magnetoelectric 
term. We will show that quantum fluctuations in such a SC im¬ 
ply an interacting topologically non-trivial phase that cannot 
be continuously deformed into the interacting topologically 
trivial one. Our claim is substantiated by a renormalization 
group analysis. This result does not follow from the classi¬ 
cal Lagrangian of the system derived earlier in Reffl5l and 
is a purely quantum effect involving the interaction between 
photons and Higgs fields. We show that due to this behavior, 
the type I regime of the TSC features a first-order phase tran¬ 
sition in the bulk and a second-order phase transition on the 
surface. This distinguishes a TSC from a topologically trivial 
superconductor in the type I regime. The latter would exhibit 
a first-order phase transition in the bulk as well as on the sur¬ 
face, provided quantum fluctuations are accounted for 2 ^ (this 
point will be discussed in detail in Section III). 


Another consequence of the topological magnetoelectric 
term is the occurrence of an anomalous Hall effect when an 
electric field is applied parallel to the surface of a TSC. Due 
to the magnetoelectric effect, a transverse current is gener¬ 
ated from a Lorentz-like force involving the relative super¬ 
fluid velocity ~ V0 and the applied electric field. The gener¬ 
ated transverse current is negative, a situation reminiscent of 
the anomalous Hall effect in superconductor^, and observed 
high-T c cuprate superconductors! 22 However, in the latter case 
the anomalous Hall effect occurs due to vortex motion induced 
by the Faraday law, and is typically a very small effect. Fur¬ 
thermore, in this case the Lorentz force acts directly on the 
vortex core, and therefore on the normal components of the 
superconductor. For this reason, it automatically leads to dis¬ 
sipation. In the case of three-dimensional TSCs, on the other 
hand, the anomalous Hall effect occurs even in absence of 
vortices and is induced solely by an external electric field via 
the magnetoelectric effect. Thus, a dissipationless anomalous 
Hall current is generated on the surface. 

The plan of the paper is as follows. In Section II we discuss 
how the Higgs mechanism works in a topological supercon¬ 
ductor. We will show that after the phases are integrated out, 
interactions between the photons automatically arise due to 
the axion term. This is in contrast to the ordinary Higg mech¬ 
anism, where the phases can be trivially integrated out by a 
gauge transformation. In Section III we discuss the role of 
quantum fluctuations and derive the effective potential on the 
surface and the renormalization group (RG) equations. We 
show that the RG equations of the topological superconductor 
cannot be connected to the ones of a topologically trivial su¬ 
perconductor. This is shown to occur as a direct consequence 
of the axion term. Finally, in Section IV we discuss the dissi¬ 
pationless variant of the anomalous Hall effect using the Lon¬ 
don limit of topological superconductors. 

II. HIGGS MECHANISM IN THREE-DIMENSIONAL 
TOPOLOGICAL SUPERCONDUCTORS 

The effective Lagrangian for a three-dimensional TSC fea¬ 
turing two Fermi surfaces is given bM 

e^O 1 

r (V - flVCTT jn F _ F F [IV 

^ eff _ 32 n 2 r w r o-T ^ r rv r 

+ Z - m2 l0;l 2 ] 

i=L,R 

- >d 2 + l^l 2 ) 2 + 2 J(<Pl<p R + <P* r (Pl), (1) 

where q = 2e is the charge of the condensate. The Greek in¬ 
dices run from 0 to 3 and F= d^A v - d v A with (Af) = 
(A 0 , A). The Lagrangian 0 corresponds to an abelian Higgs 
model with a two-component scalar field (<fc, (Pr ), where R , L 
denote ’’right” and ’’left” chiralities of the two components of 
the scalar condensate matter field, minimally coupled to the 
gauge-field (Af). In contrast to the standard Higgs model, Eq. 
Q features a so-called axion temPSl, which is the first term 
appearing in the Lagrangian above. The term is topological 
in nature and contains a scalar field (the axion) 6 = 0 L - 0 R , 
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where 0 L and 0 R are the phases of c/) L and (p R , respectively, i.e. 
(pi = \pi\ exp (iQi). In terms of electric and magnetic field com¬ 
ponents, the toplogical term reads e 2 6E • B/(4 7r 2 ), which is 
precisely the magnetoelectric form mentioned in the introduc¬ 
tory paragraphs. A Josephson coupling term oc p* L p R + P* r Pl 
accounts for the interference between the two superconduct¬ 
ing order parameter fields. This is a characteristic feature in 
superconductors with two or more components of the order 
parameter, and is absent only if prohibited by symmetries of 
the problem HUS] Below, we show that the Josephson coupling 
is generated by fluctuations, and therefore it is legitimate to in¬ 
clude such a term from the very beginning in the Lagrangian. 
Furthermore, the Josephson coupling is important for tuning 
between topologically trivial and nontrivial phases. In fact, a 
simple mean-field analysis shows that for / < 0 the Josephson 
coupling implies 0 = 0, yielding a topologically trivial super¬ 
conductor. For / > 0, on the other hand, 6 is locked to n, thus 
leading to a topologically non-trivial superconducting ground 
state. Since 6 is periodic, 6 = n corresponds to a situation 
where the time-reversal (TR) symmetry is preserved 15 . Thus, 
at the mean-field level, J = 0 separates a topologically trivial 
ground state from a non-trivial one. Thus, varying J from pos¬ 
itive to negative values induces a topological quantum phase 
transition. 

In the U( 1) Higgs mechanism the phases disappear from 
the spectrum due to spontaneous breaking of the local U( 1) 
symmetry, being transmuted into the longitudinal mode for 
the photon, which becomes gapped. Thus, only amplitude 
modes remain in the spectrum of scalar particles. The Higgs 
mechanism is equivalent to integrating out the phase degrees 
of freedom, which in the case of the Higgs model automati¬ 
cally leads to a massive gauge particle. This point of view of 
integrating out the phases is particularly appealing in the case 
where a Josephson coupling is present. However, additional 
non-linearities arise in the presence of the axion term. To see 
this, let us first consider the Higgs mechanism in Eq. 0 for 
the case where the axion term is absent. In this case we can 
simply write pj = pje l6 j / V2, j = L,R and make the shift, 


■ Ajj + 


i (p 2 L d^e L +p 2 R d^e R 


P 2 l+P 2 r 


( 2 ) 


which yields, 

£ - e 0 F _ - f F py + q — L A 2 


2 2 

PlPr 


2 (Pl+Pr ) 


(d^ef + JpLpR cos 6 


y(pi+p*)-g(pi+p*) 2 - 


( 3 ) 


Note that in absence of the axion term and for the particu¬ 
lar situation of a single scalar field component, i.e., p L = p 
and p R = 0, the above effective Lagrangian trivially reduces 
to the usual Lagrangian for the Higgs model in the unitary 
gauge. In this particular case the Lagrangian is independent 
of the phase. For the case relevant to us here, where two scalar 
fields are present, there is a term ~ ( 8^6 ) 2 remaining. Thus, 


in the absence of Josephson coupling there is still a massless 
(Goldstone) mode present in the spectrum. This occurs be¬ 
cause there are two Higgs fields and one Abelian gauge field. 
Thus, it is only possible to gauge away one phase degree of 
freedom. The phase 6 would not couple directly to the gauge 
field in the absence of the axion term. Due to the axion term, 
integrating out the phases generate direct interactions between 
photons, even if the amplitudes are assumed to be uniform. 
For the topological phase of the system occurring for J > 0, 
we have to integrate out the lowest order Gaussian phase fluc¬ 
tuations around 6 = n. In general, this renders the induced 
photon-photon interaction being non-local. If the amplitudes 
are uniform, we obtain the effective Lagrangian, 


-CHi gg s|y>o - ~~F + 

1/ e 2 \ 2 ( pI+P 2 r 

PlP 2 r 


1 r 2 , r(PL+pR> A 2 


A z + — p iva ~ r F F 
32 7C ^ aT 


2 \ 1 6n 2 J 


: J d A xy(x-x')e, Y ^ s F»\x)F*e(x)F a p(x')F y5 (x') 


- JPLPR - : f -(pI+P 2 r)~ 1(PI + Pr) 2 ’ 


( 4 ) 


where 


V(x) = 


I 


d 4 p e ip ‘ x 
(27r) 4 p 2 + m 2 ’ 


( 5 ) 


with m 2 e = J(pl/Pr + Pr/Pl )• In the long wavelength regime 
the induced photon interaction is strongly screened by the ax¬ 
ion. Hence, we have V(x - x ') « mp 2 d 4 (x - x'). The resulting 
photon-photon interaction simplifies and we obtain 27 


6 2 (P 2 l +Pr) a2 


-^Higgsl/>0 
2^ JpLpR \ 167T 2 


4 

,2 ^2 


-F + 


A z + ^—f pv(TT F F 
32;r ^ ^ 


det (/>) 


JpLpR ^ (Pl + Pr) g (Pl + 


( 6 ) 


where det^y) = (E • B) 2 . 


III. QUANTUM FLUCTUATIONS 
A. Important vanishing of a Feynman diagram 

We now turn to a crucial aspect of the topological phase 
with respect to the surface states. It turns out that in the pres¬ 
ence of quantum fluctuations, the topological surface states 
cannot be continuously deformed into topologically trivial 
ones in the long wavelength limit when crossing the critical 
point. To see, this let us assume that 6 is uniform on the sur¬ 
faces, with 6 = 7 : for the TR invariant case. Note that 6 R 
and Or are still allowed to fluctuate, with 0 L = 0/2 + S0 L and 
0 R = -6/2 + S6 r . Since e py(TT F aT = 2dP^ yXp A y F xp ), each 
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surface contains a Chem-Simons (CS) term.® Thus, assum¬ 
ing two surfaces perpendicular to the z-axis, we find that the 
imaginary time photon propagator at any surface is given by 


a Up) 


p 2 + m A 
A(p 2 ) 


PuPv , M e _ ^ 

2 ~ 7 2 e M v ^P 

p 2 p 2 +m 2 A 


(7) 


where m 2 = q 2 (p 2 + p 2 ), Mg = e 2 6>/(8 n 2 ), A(p 2 ) = (p 2 + 
m A ) 2 + MgP 2 , and the ± sign is chosen depending on which 
surface one is referring to. In Eq. 0 the transverse gauge has 
been fixed. At the phase transition to the normal state where 
m 2 A —> 0, the Feynman diagram shown in Fig. [l| which is as¬ 
sociated with Higgs scattering mediated by photons, behaves 
very differently at long wavelengths (\p\ —> 0), depending on 
whether 6 = 0 (topologically trivial) or 6 ± 0 (topologically 
non-trivial). Namely, for all 0 ± 0 we have, 

I™, lim f ^ A± ,,(/) + q)Ay (q) = 0, (8) 

bh0 m 2^ 0 J (2 ny * * 


de Gennes picture. This occurs because the photon is topo¬ 
logically gapped, despite the vanishing of the Meissner gap 
(niA = 0). Note that this result is due to the topological char¬ 
acter of the axion term, and not due to a symmetry protection. 
Indeed, since the diagram of Fig. [ljvanishes for any 6 ± 0, TR 
invariance is not required. Thus, in this context the topolog¬ 
ically non-trivial phase simply corresponds to the case where 
the axion term is nonzero. 

A continuous deformation to the topologically trivial phase 
can be done in the Higgs phase, where there are no gapless 
modes. At the critical point such a continuous deformation 
is not possible. Thus, quantum critical fluctuations in this 
system will govern topologically stable universal behavior in 
physical quantities. This has important implications for crit¬ 
ical exponents and amplitude ratios. Indeed, as we will see, 
the vanishing of the diagram shown in Fig. [T] changes signif¬ 
icantly the renormalization group (RG ))3 function associated 
to the interaction vertex between scalar fields. 


B. One-loop effective action on the surface 


whereas the result is divergent for 6 = 0. Thus, after other one- 
loop scattering amplitudes are taken into account to obtain 
the full four-Higgs vertex, we see that at the critical point the 
topological field theory cannot be continuously deformed into 
a topologically trivial one. This statement holds trivially for 
topological Bogoliubov-de Gennes superconductors. Here, 
we have shown that it also holds in the presence of quantum 
fluctuations in an interacting system, beyond the Bogoliubov- 


As is well known, the one-loop effective action is more eas¬ 
ily obtained by integrating out the quadratic fluctuations of the 
scalar fields and gauge fields. 29 We will assume that the mag¬ 
nitudes of both (p L and (pR have the same expectation value in 
the broken symmetry phase. Thus, we write (p L = e l6l2 <p + 0 L 
and (pR = e~ l6l2 <p + 0/?, where ip is uniform and 0* represents 
the fluctuation around {(pi). The effective action is therefore 
written in the form, 


^ l-loop 
eff 


+ 

i eff+ 2 


I J'd 3 x J"d 3 x' [<I) t (x)M(x - x'\ip) ®(V) + A ll (x)M /JV (x - x':ip)A v (x')\, 


(9) 


where 


= 2V[(m 2 - 27cos0)M 2 + u\<p\ 4 ], (10) 

with V being the (infinite) volume of three-dimensional spacetime, and O ' = [<p* <p* R (j) L 4>g]. The matrices M(x - x';ip) and 
M yUy (v - x'\ <p) are given in momentum space by, 


M{p, (fi) 


(p 2 +m 2 + 3u\<p\ 2 )I + (wMV^ 0 - 2J)cr x u\<p\ 2 (e Wcr * + cr x ) 

u\(p\ 2 (e~ ,e>Tz + cr x ) (p 2 + m 2 + 3m|^>| 2 )/ + (u\(p\ 2 e~ ,<Tze - 2 J)cr x 


(ID 


where / is a 2 x 2 identity matrix, while cr x and cr z are Pauli matrices, and 


Mftvip, tp) = ( p 2 + 4q 2 \ip\ 2 )S IXY - p^py - Mge^xPi- (12) 

From Eq. we see that a correction to the Josephson cou¬ 
pling has been generated by fluctuations. It has the form, 
u\(p\ 2 {e l6 (p* L <pR + e~ l6 ^ R ^> L ). This term is generated even if 
J = 0, which means that including a Josephson coupling in 
the Lagrangian is a physically reasonable assumption. This 
should be expected, since fluctuations will necessarily lead to 
an overlap between the two complex field components. Note 
that this result is valid in general for any two-component su¬ 


perconductor and is not restricted to the topological one being 
considered here. 

Integrating out the fluctuations in Eq. ([9]), we obtain, 

e ~ vu w> = J (13) 

where V is the three-dimensional spacetime volume and 
U( M, 6) is the effective potential given by, 

U(M,0) = 2 . [Tr In M^ v + Tr In m] . (14) 
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(9 = 0: 


lim 

IpI-^o 



= oc 


<9^0: 



FIG. 1. Difference in behavior for surface photon-mediated Higgs 
scattering at the critical point in the long wavelength limit. For the 
topologically trivial case (6 = 0) the corresponding Feynman dia¬ 
gram diverges. On the other hand, for the topologically non-trivial 
case (0 ± 0), the same diagram vanishes. This shows beyond the non¬ 
interacting regime that a topological superconductor cannot be con¬ 
tinuously deformed into a topologically non-trivial one in the long 
wavelength limit. 


The tracelogs can be written in more explicit form, 

U(\(p\,0) = 2[(m 2 - 2J cos 6)\(p\ 2 + u\(ff] 

+ In [p 2 + M\ a { M, 0)] + In \p 2 + M 2 2(T (M, 0)]}, (15) 

where, 

M 2 ± (\<p\, 0) = 2q 2 \<p\ 2 + ^ ± 1^1 ^M 2 + M 2 , (16) 

M 2 ± (\(p\, 0) = ot 2 + 2u\ip\ 2 ± 2 J, (17) 

M 2 ± (\<p\, 0) = m 2 + Au\(p\ 2 
± 2 ^/ 2 + m 2 M 4 - 2 Jit cos 0M 2 . (18) 

From this mass spectrum underlying the effective potential, 
we recognize p 2 + M 2 + and p 2 + M 2 _ as the would-be Gold- 
stone modes from the regime J - 0, corresponding to the ab¬ 
sence of the Josephson coupling. Note that the corresponding 
uncharged system in absence of Josephson coupling features 
three Goldstone modes. This is to be expected, since the sys¬ 
tem in this case would be 0(4) invariant. As usual, in the 
charged system there are only gapped modes, even in absence 
of Josephson coupling, as required by the Higgs mechanism. 

The momentum integral in Eq. can be easily evaluated 
using an ultraviolet cutoff A, 

U(\<p\, 0) = Uo + 2 [(m 2 - 27cos 0) M 2 + u\<p\ 4 ] 

- yE J] K(M,0)| 3 + \MxMGf 

cr=± 

+ |M 2<r (M,0)| 3 ], (19) 


where Uo is a field-independent (vacuum) contribution. We 
have absorbed a term A(3 u + q 2 )/(2n 2 ) in m 2 , since this con¬ 
tribution just corresponds to tadpole diagrams. Note that 
this simple renormalization does not actually affect the other 
terms, since the error committed by doing this is of higher 
order than one-loop. From the effective potential above we 
see that the diagram of Fig. [T] can only arise from the sum 
Zcr=+ with M 2 (\(p\, 6) given in Eq. (16). An ex¬ 

pansion in powers of \(p\ for 6 ^ 0 clearly shows that the dia¬ 
gram shown in the figure indeed vanishes when 6 ± 0. How¬ 
ever, the contributions q 2n \(p\ 2n for n > 4 become singular for 
all 6. We emphasize that this is only true when the topological 
magnetoelectric term is nonzero. On the other hand, if 6 = 0, 
we obtain, 


F \Mo.(\ip\,G = 0)| 3 = 8g 3 M 3 , (20) 

<T=± 

meaning that the effective potential is non-analytic as a func¬ 
tion of g 2 M 2 . Albeit simple, this is actually a non-perturbative 
result, since it cannot be obtained as a power series in q 2 . The 
absence of a power series involving terms g 2 "|^| 2n reflects the 
divergence of the diagram shown in Fig[l] when 6 = 0 for 
vanishing external momenta. In this case a meaningful pertur¬ 
bative evaluation of the vertex function can only be done for 
nonzero external momenta, a fact related to the infrared diver¬ 
gence arising from a massless photon. Note that in contrast to 
the 6 ± 0 case, the expansion is singular in the infrared for all 
n > 2, rather than for n > 4. However, the singularities in the 
case of a topologically massive photon is not a problem, since 
they correspond to interactions that are irrelevant in an RG 
sense. For instance, this type of infrared singularity in higher 
order vertices would also occur in a simple <£ 4 Fandau theory. 
Having 6 ± 0 turns the photon topologically massive with¬ 
out spoiling gauge invariance PD Thus, we can interpret the 
cubic contribution arising in the limit of vanishing 6 as a con¬ 
sequence of resumming up all the one-loop infrared divergent 
diagrams containing only internal photon lines. This leads to a 
non-analytic contribution to the effective potential. As is well- 
known, summing up these contributions in d = 3 + 1 in the 
case of standard scalar electrodynamics leads to a logarithmic 
term ~ g 4 M 4 In \y\ 2 in the effective potential, 30 while the cu¬ 
bic term has been obtained in the context of Ginzburg-Fandau 
superconductors for d = 3 + 0 (i.e., scalar electrodynamics in 
d = 2+1 and imaginary time). 20 In both cases these one-loop 
photon contributions lead to a fluctuation-induced first-order 
phase transition. However, it was later shown that this result 
is valid only in the type I regime, while in the type II regime a 
second-order phase transition in the so-called inverted 3D XY 
universality class arisesP 1 

In view of the above discussion, it is of interest to inves¬ 
tigate the character of the superconducting phase transition 
on the surface of a topological superconductor. The theory 
features two ingredients that are not present in the previous 
analysis of fluctuation-induced first-order phase transitions by 
Halperin et al 20 These are the Josephson coupling between 
the scalar field components, and the CS term. The case with 
Josephson coupling and in absence of CS term has been ex¬ 
amined in the Fondon limit (i.e., in the strong type II regime) 
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in Ref. j26j For this case it has been shown by means of exact 
duality arguments that a two-component superconductor with 
Josephson coupling exhibits a phase transition in the 3D XY 
universality classP 6 On the other hand, when the Josephson 
coupling is absent, it has been shown using renormalization 
group (RG) methods in Ref. 132 j that for a large enough CS 
coupling (i.e., M e in our notation), the first-order phase tran¬ 
sition is turned into a second-order one. Regarding the result 
in the London limit without CS term, we note that our anal¬ 
ysis in this Section is being done in the type I regime, since 
amplitude fluctuations play an important role in the calcula¬ 
tions above. Indeed, the theory with 6 = 0 yields a fluctuation 
induced first-order phase transition. 


C. Renormalization group analysis 


are obtained in a way similar as in Ref. [20] except that we 
use the propagator 0 with niA = 0 in the Feynman diagrams 
involving photon lines. The result is, 


dq 2 

~df 


Nq 2 
247T 


(23) 


du 

~dl 


1 + 


4 q 2 


(M + 4)., 

U H--- U , 

8tt 


where we now have, 


(24) 


6 = 



(25) 


It is a well-known fact that the CS term does not 
renormalize. 33 This is a consequence of the topological nature 
of the CS term. Indeed, since it is independent of the metric, it 
does not change under scale transformations. Thus, using this 
result and the invariance of the effective action under renor¬ 
malization, we obtain, 

Mee^A^A* = Mg^ v Xd v A^ ( 21 ) 

where A^ r = Z~ i/1 A I ‘ is the renormalized gauge field with the 
corresponding wavefunction renormalization, Za, and M ej is 
the renormalized topological mass. From the above equation 
it follows that M e , r = ZaMq. Since gauge invariance implies 
the renormalization q 2 = Z A g 2 ,®- 30 which is easily obtained 
from the vacuum polarization, it follows that 6 is a renormal¬ 
ization group invariant, 


where / = ln(m r /A) is a logarithmic renormalization scale de¬ 
fined in terms of the renormalized mass, m r , yielding the in¬ 
verse correlation length. Eq. ( [22] ) is an important result, since 
it allows one to study the critical topological behavior with 
a vanishing renormalized Josephson coupling / r , while still 
having 6 ± 0, corresponding to the critical point of the topo¬ 
logical phase transition. 

Combining both CS terms stemming from the two surfaces, 
we obtain that the axion term in the bulk does not get renor¬ 
malized either, so the result ( [22] ) also holds in the bulk. An 
interesting consequence of this analysis is that the RG flow 
of the bulk theory does not differ significantly from the one 
where the axion term is absent, which is just given by the well- 
known analysis of Coleman and Weinberg.® In this case, the 
phase transition is known to be of first-order, irrespective of 
the superconducting regime being of type I or type II. How¬ 
ever, this is not the case for the phase transition on the sur¬ 
faces. Indeed, the vanishing of Feynman graph of Fig. [ljwhen 
6 ± 0 implies that a contribution ~ q 4 is absent in the RG 
function of u. Here are q = q r lm r and u = u r lm r renormalized 
dimensionless couplings defined on the surface of a TSC. The 
one-loop RG functions for a superconductor in 2+1 dimen¬ 
sions with a CS term, and A complex order parameter fields, 


where C\j are the Chem numbers associated to the helicity 
of the A Fermi surfaces involved. 15 Note that we have not 
expanded 1/(1 + q 2 \6\/87r 2 ) in powers of q, because 6 can also 
be very large, so that the product q 2 \0\ is not necessarily small. 

It is easily seen that the above RG equations have an in¬ 
frared stable fixed point for all A, in contrast to the analysis by 
Halperin et al. for the non-topological superconductor, where 
infrared fixed points are only found for A > 183.® However, 
there is a stability condition involving 6 that has to be fulfilled 
in the case of a TSC. It is obtained by considering the critical 
correlation function of the superconducting order field com¬ 
ponents, 

<0 i( *)0}(O»~ A (26) 

where at one-loop, 



which implies the inequality q > -1. The latter inequality is 
fulfilled provided, 

|0|> |(4VV-A0, (28) 

and we see that for A > 16 a quantum critical point is obtained 
for all values 6 , showing that at least sixteen Weyl fermions 
are necessary to have a quantum critical point. Values of 6 vi¬ 
olating the inequality q > -1 correspond to a situation where 
a continuum limit cannot be defined and is therefore unphys¬ 
ical. Thus, in order to have a physically meaningful phase 
transition, the inequality ( [28] ) has to be satisfied. We note that 
the lower bound for \6\ is larger than n when A = 2. Therefore, 
the TR symmetric value 6 = n obtained at the mean-field level 
does not produce a second-order phase transition when quan¬ 
tum fluctuations are accounted for. The one-loop RG predicts 
that a second-order phase transition occurs for 6 = n only if 
A >10. 

It is tempting to relate the critical value of A to the Z \6 
classification. 19 However, at this stage it would be too specu¬ 
lative, since higher order results may affect the values of A for 
which critical points obey the inequality q > -1. 
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As a final remark on the RG analysis, let us comment on 
another possible renormalization scheme allowing to contin¬ 
uously connect the cases 6 = 0 and 6 ± 0. The infrared di¬ 
vergences stemming from the photon propagator for 6 = 0 
requires defining the renormalized four-point vertex at exter¬ 
nal nonzero momenta. This yields a renormalization scale yu 
that replaces m r in the RG flow. Since the diagram of Fig. 
[T] does not vanish for nonzero external momenta, it turns out 
that a q 4 term would be generated in the RG f3 function of u , 
even when 6 ± 0. This q 4 term would have a ^-dependent co¬ 
efficient allowing to smoothly connect the result to the known 
RG equations of a topologically trivial superconductor in the 
limit 6 —> 0. We would find once more that a second-order 
phase transition occurs for large enough 6 and a stability cri¬ 
terion would follow from the inequality q > - 1 P 2 The latter 
inequality would yield in this case values of 6 leading to a 
negative coefficient of the q 4 term, with values of 6 yielding 
a positive coefficient of q 4 violating the condition q > - 1. 
This implies that the RG f3 function of u can actually not be 
continuously connected to the 6 = 0 regime, since to this end 
it would be necessary to enter a regime where the continuum 
limit is not even defined. 


D. Summary of the phase structure 

There are two important consequences of quantum fluctua¬ 
tions as unveiled by the analysis in this Section. First, we note 
that it is not possible to reach the topologically trivial phase 
from the topologically nontrivial one within the RG. Simply 
taking the limit 6 —> 0 does not recover the RG flow of topo¬ 
logically trivial superconductors, while this limit can be real¬ 
ized classically. Second, and most importantly, the phase tran¬ 
sition in the bulk is always a first-order one, while a second- 
order phase transition is possible on the surface. This is not 
the case for the topologically trivial superconductor, where a 
second-order phase transition on the surface is only obtained 
for sufficiently large N. 


IV. VORTEX-FREE ANOMALOUS HALL EFFECT 


We next turn to the Meissner effect aspects of a TSC, which 
as we will now show, implies an anomalous Hall effect even 
in the absence of vortices. This is more conveniently done by 
rewriting the Lagrangian in a London limit exhibiting explic¬ 
itly electric and magnetic fields, i.e., 

i, = i(E 2 -B ! ) + 0EB 

+ \ y 2 - qAf,) 2 + Jp L p R cos 0 

^ i=L,R 

- y(pi + p«)-^(pi+p«) 2 - (29) 

From the effective Lagrangian we obtain that the electric 
displacement and magnetic fields are given respectively by 


D = E + e 2 6B/n and H = B - e 2 6E/n , while the supercon¬ 
ducting current is given by, 

j = q\pl(ye L - qA) + p 2 R (ye R - 9 a)]. ( 30 ) 

From Eq. ( [30] ) we obtain the usual London equation in ab¬ 
sence of vortices, V x j = -(l//l 2 )B, where A 2 = l/m is the 
square of the penetration depth. Thus, the Maxwell equation 
in the presence of the axion field, 

e 2 

V x B = j + d t E + — (ye x E + <9,0 B), (31) 

71 

yields the equation determining the London electrodynamics 
of the TSC in the form, 

e 2 

d 2 B - V 2 B + m \B = — [V x (V6 x E) 

71 

+ V x (d t 6 B)]. (32) 

For the axion field we obtain the equation of motion, 

d 2 0- V 2 6> + m 2 sin6> = + -Ue-B. (33) 

4n ' \P L P R ) 

In the low frequency regime and in absence of vortices, the 
London equation ( [32] ) simplifies to, 

e 2 

- V 2 B + ml B = — [V x (V6» x E)], (34) 

71 


while the current satisfies, 


e yn 

- V 2 j + m\j = -—(V0 x E). (35) 

71 

The London equation for the electric field is unaffected by the 
axion term, retaining its traditional form, -V 2 E + m^E = 0. 
This result is closely related to the fact that the electromag¬ 
netic energy density does not contain a magnetolectric term. 

We now consider a solution with a simple geometry, 
namely, a semi-infinite TSC (z > 0) with a surface at z = 0 at 
an external electric field E 0 = E 0 x parallel to the surface. We 
obtain, 


d 2 e 2 • 

—r + m G sin 6 : 

dz z 


Hk* k' EAz,BM - 


where E x (z) = E$e mAZ , and 

d 2 B 


e 2 d 


+ m 2 A B x = - — j- 
dz 2 n dz 


' , d0 

E x(Z) — 

dz 


(36) 


(37) 


d 2 jy 


+ m AJy 


e -^AE x {z)^-. 

dz 


dz 2 n 

The solution for Eq. ( |38l > in terms of the axion is, 

^2^,2 ; 


jy(z) 


e 2 m 2 A E 0 


m t 

2 m A 


/~*oo 

I *■ 


(38) 


-2 m A z' 


0(z) 


( 39 ) 
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FIG. 2. Induced anomalous Hall current j y (z) for m e /m A = 
1,1/2,1/4,1/6. 


where 0(0) = n. Since E • (V x j) = -A 2 E • B, we obtain the 
following relation, 


e 2 m 2 A En 


-B x (z) = m A j y (z) + — ^e~ mAZ [0(z) - 0( 0)], 


(40) 


which implies m A B x (0) = j y ( 0). Thus, we find that the 
usual boundary condition of the London theory, d j y /dz\ z =o = 


m A j y (P), is obviously fulfilled by the solution (39) in the pres¬ 
ence of the axion field. However, the Maxwell equation pT) 
in the static regime implies a boundary condition that deviates 
from the standard one in the London theory of non-topological 
superconductors, 


dBx 

dz 


= ^A^x(O) + 

z=0 


e 2 Eo d6 
n dz 


z =o 


(41) 


From Eqs. ( [39] ) and ( [40] ) we see that an approximate solu¬ 
tion can be obtained by considering terms proportional to e A 
as being of higher order, which amounts to approximating Eq. 
( [36] ) as being homogeneous. In this case we can use the do¬ 
main wall solution 0(z) = n + 2 arcsin[tanh (m e z)] in Eq. ( [39] ), 
which yields j y (z) explicitly. The explicit solution for j y (z) 
with m e ± m A in terms of hypergeometric and Lerch transce- 
dents, is not very illuminating. Instead, we plot it in Fig. [2] for 
four different values of the ratio m e lm A . It has a negative sign, 
just like in the case of the anomalous Hall effect in high-T c 
superconductor! 22 . As emphasized in the introductory para¬ 
graphs, the anomalous Hall effect in non-topological super¬ 
conductors has a quite different origin from the one discussed 
here. In three-dimensional TSCs the anomalous Hall current 
arises independently of vortex motion and is associated with 
a dissipationless current. 

For me = m A (blue curve in Fig. [2]), the expression for j y (z) 
does not involve special functions, reading, 


jy(z)=—^{(n/2)e~™-2 
7T 

- 2[e~ mAZ arctan(e mAZ ) - e mAZ arctan(e _mAZ )]}. (42) 


In order to make connection with the quantization of Hall 
conductivity in the normal state, it is instructive to integrate 
the current density over z e [0, oo) to obtain the surface current 
density, 


•surf 

Jy 


Jo 


dzj y (z) = —E 0 0(0) 


e 2 E 0 m A r 

x Jo 


dz(e 


~ m AZ _ -2m a z 


)0(z). (43) 


We obtain, 


•surf _ 
Jy ~ 


e 2 E« 


2 n 


m e + m A \ I m g + 2 m A 


+ *( 


3 me + 2 m A 
4 me 


■ 2 ^ 


4 me 
3 me + m A 
4 m e 


(44) 


where if/(z) = T\z)/T(z) is the digamma function. The nor¬ 
mal state corresponds to m A lm e —> 0, which yields 7 ^ urf = 
~(e 2 /2)Eo , i. e., the half-quantum of the quantized Hall con¬ 
ductivity. 


V. CONCLUSION 

In conclusion, we have shown that due to quantum elec¬ 
tromagnetic fluctuations, the Higgs mechanism in three- 
dimensional TSCs implies a robust topological state of mat¬ 
ter, since its RG flow cannot be continuously deformed into 
the RG flow of a topologically trivial one. This is an exam¬ 
ple of a topological state that is protected due to the coupling 
of phase and electromagnetic fluctuations via the axion term, 
with TR symmetry not being required. In fact, TR can be 
spontaneously broken by quantum fluctuations. In this con¬ 
text, we have also shown that a second-order quantum phase 
transition happens on the surface of a TSC, while its bulk un¬ 
dergoes a first-order phase transition. Without the axion term 
a first-order phase transition would happen both in the bulk 
and on the surface, provided the superconductor is in the type 
I regime. 

Another aspect of the Higgs mechanism we have studied 
is the influence of the axion term in the Meissner effect. We 
have found that the gradient of the axion field on the surface 
induces a transverse supercurrent. In the low frequency limit 
this implies a London regime leading to the generation of an 
anomalous Hall current with a negative sign. This anoma¬ 
lous Hall current is dissipationless and is the consequence of 
a Lorentz-like force involving the relative superfluid velocity 
which is simply given by the gradient of the phase difference 
between the chiral superconducting components. 
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